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The equilibrium sliort-range order (SRO) in Cu-Pd alloys is studied theoretically. The evolution 
of the Fermi surface-related splitting of the (110) diffuse intensity peak with changing temperature 
is examined. The results are compared with experimental observations for electron-irradiated sam- 
ples in a steady state, for which the temperature dependence of the splitting was previously found 
in the composition range from 20 to 28 at.% Pd. The equilibrium state is studied by analysing 
available experimental and theoretical results and using a recently proposed alpha-expansion the- 
ory of SRO which is able to describe the temperature-dependent splitting. It is found that the 
electronic-structure calculations in the framework of the Korringa-Kohn-Rostoker coherent poten- 
tial approximation overestimate the experimental peak splitting. This discrepancy is attributed to 
the shift of the intensity peaks with respect to the positions of the corresponding reciprocal-space 
minima of the effective interatomic interaction towards the (110) and equivalent positions. Com- 
bined with an assumption about monotonicity of the temperature behaviour of the splitting, such 
shift implies an increase of the splitting with increasing temperature for all compositions considered 
in this study. The alpha-expansion calculations seem to confirm this conclusion. 
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I. INTRODUCTION 

Almost a decade ago, Kulik et al% published their ex- 
perimental results on electron diffraction from irradiated 
Cu-Pd alloys. In that study samples with 20, 22, 24 and 
28 at.% Pd were maintained by high-energy electron ir- 
radiation in a steady disordered state away from their 
thermal equilibrium state at temperatures between 200 
and 400 K. At these temperatures and compositions equi- 
librium Cu-Pd alloys exhibit long-range order; the disor- 
dered state occurs only at much higher temperatures. In 
the equilibrium disordered state the intensity of diffuse 
scattering from Cu-Pd alloys with more than about 15 
at.% Pd is characterised by the fourfold splitting of in- 
tensity peaks located at the (100)™ tilO) and equivalent 
positions in the reciprocal space.ErLl The resulting dif- 
fuse intensity distribution has maxima at the (IqO) and 
equivalent positions (Fig. |l|); the value of q increases with 
increasing Pd concentration. This fine structure of dif- 
fuse scattering is caused by the atomic short-range or- 
der (SRO) and is a result of the indirect interaction of 
alloy atoms through conduction electrons in a situation 
when an alloy has reasoriahly well-defined Fermi surface 
with relatively flat areas .I'EI In this case the correspond- 
ing minima of the effective pair interatomic interaction 
in the reciprocal space are also split, and their separation 
is related to the wavevector 2ki? which spans these flat 
areas of the Fermi surface. 

Similar splitting of the diffuse intensity peaks was ob- 
served in the nopcquilibrium steady disordered state un- 
der irradiation.Q In addition to the expected concentra- 
tion dependence of q, its variation with irradiation tem- 
perature was found. Even more curious was the qualita- 
tive change of the temperature dependence of the split- 
ting with concentration: q decreased with increasing tem- 
perature in the case of 20, 22 and 24 at.% Pd, but this 



trend was reversed for the alloy with 28 at.% Pd for 
which an increase of the peak separation with temper- 
ature was found. At the same time, there was an in- 
crease of the scattering intensity with increasing temper- 
ature for all four compositions, contrary to the case of 
alloys at equilibrium where the intensity decreases with 
increasing temperature. A qualitative explanation of the 
temperature dependence of the splitting was proposed as 
follows. Firstly, the behaviour of the equilibrium SRO 
diffuse intensity in the case of the exactly solvable one- 
dimensional Ising model with competing antiferromag- 
netic nearest- and next-nearest-neighbour interactions 
was studied. It turned out that the peak positions var- 
ied with temperature. This result is in contrast with the 
mean-ifield-related Krivoglaz-Clapp-Moss (KCM) treat- 
mentoM which predicts temperature-independent peak 
positions at the minima of the interaction. As the inten- 
sity increased with decreasing temperature, the peak po- 
sitions shifted towards the wavevector of the correspond- 
ing ground state. A similar result was obtained earlier 
for the two-dimensional ANNNI model using the cluster 
variation method.EJ It was concluded that the tempera- 
ture dependence of the peak positions is a phenomenon 
which cannot be understood in the framework of mean- 
field theory. Secondly, the assumption was made that 
the behaviour of the diffuse intensity in irradiated Cu- 
Pd alloys was analogous to that of the equilibrium one- 
dimensional model. The only qualitative difference be- 
tween the two cases was the opposite roles played by tem- 
perature. Based on this "inverse temperature hypothe- 
sis" , the conclusion was drawn that one might expect to 
find the increase in q with temperature for the Cu-Pd 
alloy system at equilibrium. Here it may be added that 
the reversal of the temperature behaviour of the peak 
splitting could be expected according to this hypothesis 
as concentration increases, from the increase with tem- 
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FIG. 1. Schematic reciprocal-space picture of scattering 
from disordered Cu-Pd alloys. Large dots represent the Bragg 
reflections. Small dots correspond to the split diffuse intensity 
peaks. 



splitting with changing temperature was found. 

The idea of the present study is to gain further insight 
into the behaviour of Cu-Pd alloys under irradiation by 
studying theoretically the evolution of the diffuse peak 
splitting with changing temperature in these alloys at 
equilibrium. More exactly, the aim is to find out whether 
the splitting increases with temperature in the range from 
20 to 24 at.% Pd as the inverse temperature hypothesis 
implies and whether this behaviour is reversed as con- 
centration increases to 28 at.% Pd. Starting from the 
experimental SRO diffuse intensity measured at a par- 
ticular temperature, we first solve the inverse scatter- 
ing problem and calculate the effective interaction which 
is assumed to be pairwise and temperature-independent. 
Than the direct problem is solved and the self-energy and 
diffuse intensity at different temperatures are calculated. 
The underlying theory of SRO is described in Sec. ||. 
Sec. Ill considers data for Cu-Pd alloys existing in the 



literature. Finally, the results are discussed in Sec. 



II. ALPHA-EXPANSION THEORY OF SRO 



peraturc for 20, 22 and 24 at.% Pd to the decrease for 28 
at.% Pd. 

Last year, such an increase of the peak separation 
with tempfttature was observed at equilibrium by Re- 
ichert et-ai£3 for the disordered CusAu alloy. Moss and 
Reichertll3^t3 found the same behaviour by anaU^sing the 
Monte Carlo simulation results of Roelofs et alS3 for the 
Cu-I4.4 at.% Al alloy. In the latter work the inverse 
Monte Carlo pair interactions were determined from the 
experimental diffuse intensity at a single temperature and 
subsequently used to generate the Monte Carlo intensities 
at other temperatures. The theory of the temperature dp-, 
pendence of the splitting was proposed by Tsatskisjlj'EEl 
it identifies the wavevector dependence of the self-energy 
of the pair correlation function (PCF) as the origin of 
this effect. The self-energy S(k) and the interaction term 
2/3y(k) enter the expression for the SRO diffuse intensity 
on an equal footing (Eq. ( p. la ) below). In the KCM ap- 
proximation the fact that the self-energy is a function of 
k is ignored. Apart from the observed increase in q with 
temperature, the possibility of the opposite behaviour, 
i.e., the decrease of the peak separation as temperature 
increases, was predicted. This is exactly what should 
be expected under equilibrium conditions for the Cu-28 
at.% Pd alloy, if the inverse temperature hypothesis is 
valid, although such temperature dependence was never 
seen experimentally. The possibility of the reversal of 
the temperature dependence of the peak splitting with 
increasing concentration seems to be indicated also by 
thoj-esults of the Monte Carlo simulations of Ozoliiis et 
al.cJ for the first-principles alloy Hamiltonian with pair 
and multiatom interactions (25 at.% Pd) aiid the X-ray 
scattering measurements of Reichert et alE (29.8 at.% 
Pd). In both cases no (or a very small) change of the 



We start by describing the theory of SRO which leads 
to the t emp erature-dependent peak splitting and is used 
in Sec. |^ to relate the SRO diffuse intensities at dif- 
ferent temperatures. This theory isJaased on the alpha- 
expansion (AE) for the self-energy;Ej the self-energy is 
the only unknown quantity in the otherwise formally ex- 
act expression for the SRO intensity. The AE is the ex- 
pansion in powers of SRO parameters aimm hence the 
name. It was proposed as a generalisation of another ap- 
proach to the calculation Oit-tbei self-energy, the gamma- 
expansion method (GEM),t3~H3 to deal with distant in- 
teractions which are essential in the case of the Fermi 
surface-related splitting. The complete set of the AE 
equations has the form 
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(2.1a) 



c(l - c) [-I](k) -f 2/3y(k)] ' 
S(k) = Eooo + J2 ^ 

(k) , (2.1b) 

Imn^OOO 

^Imn = aaLn + ^"fmn i ¥= 000 , (2.1c) 

aooo = ^ y dklsRoi^) = 1 , (2. Id) 

aimn = ^ y dk /si^o (k) A;„„i (k) . (2.1e) 



In Eqs. (2T) k is the wavevector, IsRoi^) is the SRO 
diffuse intensity in Laue units, c is the concentration, 
S(k) is t he s elf-energy of the PCF G (the latter is defined 
by Eq. (|2j) below), /3 = l/T, T is the temperature in 
energy units, and V^(k) is the Fourier transform of the 
pair ordering potential 
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The potential V^^ corresponds to the interaction be- 
tween an atom of type a at site i and an atom of type /3 
at site j. Further, aimn, ^imn-, Zimn and 



Aim„(k) = Z;J„ ^ exp(zkr) 



(2.3) 



are the SRO parameter, matrix element of the self-energy, 
coordination number and shell function for the coordina- 
tion shell Imn, respectively, while aooo and Sqoo are the 
corresp ondin g diagonal matrix elements. The summation 
in Eq. (2.1b) is performed over all coordination shells, 
whereas that in Eq. (|2.3| ) is over the lattice vectors r be- 
longing to th e coordination shell Imn. The integration 
in Eqs. ( 2. Id ) and ( 2.1e ) is carried out over the Brill ouin 
zone of volume fi. Coefficients a and b in Eq. (2.1c) are 
functions of concentration, 



b = 



(l-2c)^ 
2[c(l-c)]2 ' 

[1 - 6c(l - c)]' - 3(1 - 2c)* 
6[c(l-c)]3 



(2.4a) 
(2.4b) 



The SRO parameters a are proportional to the corre- 
sponding matrix elements of the PCF G, 

Gf^ = Gf^ = ~G^^ = c(l - c)a., , (2.5) 

the definition of the PCF being 



g: 



where is the occupation number, 



1, atom of type a at lattice site i 
'0, otherwise , 



(2.6) 



(2.7) 



and angular brackets denote statistical averaging. 

The meaning of Eq s. (^ is as follows. The first of 
Eqs. ( p.ld| ) and Eq. (2.1e) are the consequences of the 



fact that aij is the back Fourier transform of the intensity 



IsRoi}^)- Eq. (2. It) is the relation between the direct- 
and r ecipr ocal-sp ace representations of the self-energy. 
Eqs. ( pTb| ) and ( pT^ are written in coordination shell 
notations— The second of Eqs. (2. Id) is the well-known 
sum rulecll which reflects the property 



(2.8) 



of th e oc cupation numbers follo wing from their defini- 
tion (2/7). The less obvious Eq. (bJa ) is one of the pos- 
sible forms of the Dyson equatiorEJwhich is satisfied by 
the PCF (b.fj) i-t|his issue is discussed in considerable de- 



tail elsewhere m The key equation is Eq. ( |2.1c| ) which 
closes the set of Eqs. (2.1) by expressing the off-diagonal 
part of the self-energy in terms of the SRO parameters. 
Its right-hand side is, in fact, two first non-zero terms 
of a series expansion of Tjimn in powers of the SRO pa- 
rameters. The latter are almost always sufficiently small. 



which justifies thcj^exnansion. These two terms were pre- 
viously calculatedEj'Ey in the framework of the GEM us- 
ing self-consistent renormalization of the bare propaga- 
tor [[3V)~r-^ in the generating functional for correlation 
functions J13 The resulting expansion for the matrix el- 
ements of the self-energy was in powers of the matrix 
elements of the fully dressed propagator. This propaga- 
tor is the PCF (2^), and its matrix elements are there- 
fore proportional to the corresponding SRO parameters. 
Thus, Eqs. (|]^) form the set of self-consistent equations 
for the matrix elements of the self-energy (alternatively, 
Sooo and aimn, Imn ^ 000, can be used as independent 
variables) and constitute the closed-form approximation 
for SRO. A pa rticular AE approximation is defined by 
using Eq. ( 2.1c ) for only a finite number of coordination 
shells and neglecting all other matrix elements of the self- 
energy. Another sequence of the AE approximations can 
be generated in the same way by taking into account 
only the lowest-order (quadratic) term in the AE expan- 
sion for Sjmn and ignoring the third-order contribution. 
For the rest of the paper both terms (as in Eq. (2.1c)) 
will be used. The AE is expected to be at least as ac- 
curate as the GEM, and the latter was used successfully 
in dealing wi th . b |at|i direct and inverse diffuse scatter- 
ing problems ,113123^0 providing reliable results at almost 
all temperatures. The zero-order approximation of the 
AE i^he well-known spherical model (SM) for correla- 
tions ,l3 alsoinown under the name of the Onsager cavity 
field theory.E3 In the SM the self-energy is diagonal, i.e., 
wavevector- independent; the single non-zero matrix ele- 
ment Sqoo is a function of temperature and concentration 
and is determined from the sum rule (2. Id). 

In order to use Eqs. (2.1) for calculating the evolution 
of the diffuse intensity with temperature and, in particu- 
lar, the temperature dependence of the peak splitting, it 
is necessary to have information about the interaction V . 
It is assumed here that the interaction does not depend 
on temperature in the relevant temperature intervals; on 
the other hand, it is clearly concentration-dependent, so 
that a separate interaction set is needed for each alloy 
composition. We start from the set of the experimental 
SRO parameters and calculate the AE interaction in the 
reciprocaL^pace by solving the inverse diffuse-scattering 
problemO^EJ This interaction can then be used for calcu- 
lation of diffuse intensities at different temperatures and 
possibly, with much less confidence, at slightly different 
concentrations . 



To solve the inverse problem, we rewrite Eq. ( 2.1a ) as 
an expression for the interaction: 
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E(k) 



(2.9) 



The SRO diffuse intensity here is recalculated from the 
set of the experimental SRO parameters: 

/sj?o(k) = l+ ^ 

(k) . (2.10) 
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In Eq. (2.1C) the sum rule (2. Id) was used; otherwise, 
it is just the Fourier transformatio n wri tten in coordina- 
tion s hell n otatio ns, s imilar to E g. (^.Ibl ) . Substitution of 
Eqs. ( ^.lb|) and (|2J^ ) into Eq. (|2.9D shows that the only 
quantity in the resulting expression for VaeO^) which 
is needed to be expressed in terms of the SRO param- 
eters (or, equivalently, the SRO intensity) is the diago- 
nal part Sooo of the self-energy. The off-diagonal part 
of S is alre ady an explicit function of the SRO param- 
eters (Eq. ( |2.1c )). To find Eqoo, we integrate Eq. 
over the Brillouin zone; this integration gives the diag- 
onal direct-space matrix element of the integrand, as in 
Eq. (2. Id). The interaction V is an off-diagonal matrix 
in the direct space because of the absence of the self- 
interaction. Therefore, after the integration the left-hand 
side of Eq. (2.9) is zero and, as a result. 
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c(l - c) 



(2.11a) 
(2.11b) 



Thus, Eq. (2_^) for the AE interaction can be written as 



VAB(k) = FsA/(k) + ^Sod(k) , 
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FIG. 2. CoHicentration dependence of the peak splitting 
q as measuredtl in the X-my diffraction experiment (circles, 
dashed line) and estimatedcl fromr4he results of the KKR-CPA 
electronic-structure calculation£j (squares, dotted line). The 
data am taken from Table 1 in Ref. BI The original KKR-CPA 
result£3 (triangles, solid line) are also shown. Straight lines 
connecting symbols are for the eye guidance only. Note that 
all the data were originally given in units of the distance be- 
tween the (000) and (200) positions (equal to 2 r.l.u.) for the 
separation m = \f2q between the adjacent peaks. 



where Eod(k) is the Fourier transform of th e off-diagonal 
part of the self-energy defined by Eq. (2.1c), and 



VsM(k) 



T 



2c(l - c) 



VsRo{^) - (^sAo)] (2-13) 



is the interaction obtained in the framework of the SM, 
i.e., in the zero-order AE approximation in which the off- 
diagonal part of the self-energy is zero. Note that, when 
compared with the intefiae|tion resulting from the KCM 
expression for Isno^P^ 



VR-cAf(k) 



T 



2c(l - c) 



[^sAo(k)-l] 



(2.14) 



the SM interaction differs, at given c and T, only by the 
constant subtracted from the inverse intensity. There- 
fore, the off-diagonal direct-space interaciions are identi- 
cal in the KCM and SM approximati ons £j However, the 
KCM formula violates the sum rule ( 2. Id ), thus leading 
to the appearance of the unphysical self-interaction 
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2c(l 



(2.15) 



while in the SM, according to Eq. (|2.13| ), this matrix 
element is zero. Returning to Eq. ( |2.12| ), every term in 
its right-hand side is expressed at this stage in terms of 
experimental data, and VA_B(k) can be easily calculated. 



III. AVAILABLE DATA 

We now consider previously published experimental 
and theoretical results for equilibrium Cu-Pd alloys in 
the discussed range of concentrations (20 to 30 at.% Pd). 
These results are of two types. Firstly, the electron and 
X-ray diffraction data and the results of the Korringa- 
Kohn-Rostoker coherent potential approximation (KKR- 
CPA) electronic-structure calculations are available for 
the concentration dependence of the peak splitting. Sec- 
ondly, for several alloy compositions large sets of the SRO 
parameters were determined by the Fourier inversion of 
the experimental SRO diffuse scattering intensities. The 
latter type of data is used as an input for the calculations 
of the kind described in Sec. ^ 

The peak separation q was measured at equilibrium 
for yatious concentr8|tions and temperatures using elec- 
tronlicl and X-rayETIJ scattering. The splitting was ob- 
served for alloys with more than about 15 at.% Pd, and 
it increased monotonically with increasing Pd content. 
Thou^t.very good agreement was noted by Ciyorffy and 
Stock£3 (GS) between the electron-diffractiona and their 
KKR-CPA results, the calculated values of q were sys- 
tematically slightly higher than the experimental ones. 
The discrepancy became noticeably larger in more recent 
measurements. In particular, Saha et aln compared their 
X-ray results for several compositions with the estim|a=. 
tions they made from the GS KKR-CPA calculations .EZl 
They found that the experimental splitting was smaller 
in all cases. The difference in q ranged from 0.014 to 
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FIG. 3. The same as in Fig. bL but all available experi- 
mental data are presented (Ref. Tl - open squares, ca. 800 
K; Ref. ^ - open triangles, 773-893 K; Ref. ^ - open circles, 
ca. 700 K; Ref. | - plusses, 1073 K; Ref. | - crosses, 1023 K; 
Ref. P - asterisk, 773 K; Ref. |^ - open diapiond, ca. 700 K), 
m comparison with the KKR-CPA resultfO (filled triangles, 
solid line) . The estimations made in Ref. |^ are not shown. 



0.034 reciprocal lattice units (r.l.u.); 1 r.l.u. is the dis- 
tance between the (000) and (100) positions. Their find- 
ings (Table 1 in Ref. ^) are shown in Fig. ||, together 
with the original GS results for different concentrations 
which were read off Fig. 2 in Ref. Surprisingly, all 
the estimated values, presumably calculated by interpo- 
lating the results of Ref. |2^, lie below the GS line; the 
reason for this is not clear. As a result, the disagreement 
between the experimental and theoretical values of q is 
even more pronounced than it was reported in Ref. |^. 
Fig. ^ compares the KKR-CPA results with the collec- 
tion of all experimental data for the splitting known to 
the authors. It is seen that all the experimental points 
are located below the GS line. We propose an explana- 



tion for this discrepancy which is given in Sec. IV 



Sets of the SRO parameters were obtained in the con- 
sidered concentration interval for 21.8, 28.5 (Ref. ||) and 
29.8 (Ref. H) at.% Pd in X-ray experiments. The sam- 
ples were annealed at some temperature corresponding to 
the disordered phase and then quenched. Hereafter these 
alloys will be referred to according to their numbers in 
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FIG. 4. Profiles of the recalculated SRO diffuse intensities 
for the alloys 1 (a), 2 (b), and 3 (c) along the (hlO) line. 
Variable k is the component of the wavevector k = (fc, 1, 0). 
Note that there is no splitting of the (110) peak for the first 
two alloys. 



TABLE I. Data for three Cu-Pd alloys for which sets of the 
SRO parameters are available: T is the annealing tempera- 
ture, Na the number of the SRO parameters in the set, aooo 
and g'^^'' the experimental values of aooo and g, respectively, 
q^'^'^ corresponds to the recalculated intensity (see text). The 
splitting q is measured in r.l.u. 



No. 


at.% Pd Ref. 


r, K 




exp 

"ooo 


^exp 




1 


21.8 






1023 


78 


1.018 


0.059 





2 


28.5 






1023 


78 


1.014 


0.151 





3 


29.8 






773 


72 


1.786 


0.184 


0.162 



Table | which contains data used in the subsequent dis- 
cussion. The splitting of the experimental (110) intensity 
peak was detected for all three compositions. The SRO 
parameters for large number of coordination shells were 
calculated by Fourier-transforming the SRO part of the 
measured diffuse intensity after having separated it from 
other intensity contributions. We recalculated the SRO 
diffuse intensities for these three alloys using tables of the 
SRO parameters given in Refs. ^ and ^ and the theoret- 
ical value aooo = 1 instead of the experimental values. 
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FIG. 5. Recalculated SRO intensity for the alloy 2 along 
the (hOO) line showing ranges of negative values. Variable k' 
is the component of the wavevector k = {k' , 0, 0). 



Such substitution leads to a simple shift of the intensity 
and does not change its shape. Surprisingly, no split- 
ting of the (110) peak was found in the recalculated SRO 
intensities for the alloys 1 and 2 (Figs, ^(a) and ^(b)); 
in the case of the alloj2-.2 this circumstance has already 
been noted elsewhere.ES In addition, negative values of 
the recalculated SRO intensity were found for the alloy 2 
(Fig. H). The origin of all these inconsistencies is prob- 
ably the insufficient accuracy and/or number of the cal- 
culated SRO parameters. Contrary to these two cases, 
the recalculated SRO intensity for the alloy 3 shows the 
experimentally observed splitting (Fig. ^(c)). However, 
the corresponding value of q is noticeably smaller than 
the experimental result (Table |). As for the first two al- 
loys, we find that the splitting tends to decrease after the 
recalculation. The accuracy of the recalculated intensity 
seems to be better for the alloy 3 as far as the magnitude 
of the splitting is concerned, though the deviation of the 
integrated intensity aooo from unity, which often serves 
as an accuracy criterion in diffuse-scattering experiments, 
is much larger than in two other cases (Table |). 



IV. RESULTS AND DISCUSSION 

We assume that the discrepancy between the experi- 
mental and theoretical values of q discussed in Sec. [II is 



the result of the shift of the intensity peak position with 
respect to the position of the,corresponding minimum of 
the interatomic interaction.E3 In other words, quantities 
which were measured and calculated MBve not the same. 
Indeed, what GS actually calculatedEZi using the KKR- 
CPA method were the Fermi surfaces and, in particular, 
the Fermi wavevectors kj? along the (110) direction for 
different concentrations. These Fermi wavevectors were 
subsequently used to calculate the 2ki?-related separa- 




FIG. 6. Shift of the intensity peak position as a result of the 
wavevector dependence of the self-energy. Behaviour of the 
SRO intensity, self-energy and interatomic interaction along 
the (hlO) line is shown schematically. The self-energy profile 
is as found for the three Cu-Pd alloys discussed in the text 
(see Fig. |^ below). Dashed lines indicate positions (left to 
right): (110), intensity peak, minimum of the interaction. 



tion m — \p2.q between the adjacent minima of the inter- 
action T^(k). Since the mean-field (KCM) description of 
correlations was chosen, the resulting separation between 
the intensity peaks was the same. However, it is generally 
different from the separation between the l^(k) minima 
and depends on temperature because of the tetaperature- 
dependent shift of the intensity peak position.ta The shift 
itself is the consequence of the wavevector dependence of 
the self-ei rergy (Fig. This can be easily seen from ei- 
ther Eq. ( ^.la| ) (the direct problem) or Eqs. ( p^), ( |2l^ ) 
(the inverse problem). Consider, e.g., Eq. (2.1a); the 
-^Sflo(k) peak positions are determined by the condition 
V/sflo = 0, which leads to 



2VF = TVE 



(4.1) 



while the positions of the l^(k) minima are obtained from 
the equation VV^ = 0. It is clear from Eq. (4.1) that the 
extrema of /si?o(k) away from the special points are, in 
general, different from those of ^(k). On the other hand, 
if the approximate self-energy is k-independent (as in the 
KCM or the SM approximations), then the two equations 
coincide and the intensity peak is not shifted. Fig. ^ 
shows that the KKR-CPA calculations overestimate the 
experimental peak splitting everywhere in the range of 
concentrations from 20 to 30 at.% Pd. In the frame- 
work of the suggestion about the shift of the /si?o(k) 
peak being the origin of the disagreement between ex- 
periment and theory, this means that the intensity peaks 
are shifted towards the (110) position. 

Based on this assumption, it is now possible to pre- 
dict the temperature behaviour of the splitting if another, 
sufficiently reasonable assumption is made. We assume 
that the temperature dependence of the splitting is al- 
ways monotonic (a non-monotonic behaviour was never 
observed experimentally). If this assumption is correct. 
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FIG. 7. AE effective pair interactions VAE(k) for the al- 
loys 1 (a) and 3 (b) along the (hlO) line. 



then the direction of the shift at a particular tempera- 
ture value can be related to its temperature dependence. 
At high temperatures corrections to the KCM approx- 
imation are small, and the absolute value of the shift 
tends to zercw-decreasing at least as with increasing 
temperatureJl3 Therefore, in the case of the monotonic 
behaviour of the splitting the direction of the shift is the 
same at any temperature; the splitting increases with 
temperature, if the shift is towards the (110) position, 
and decreases otherwise. For Cu-Pd alloys this would 
mean that the splitting increases with temperature for 
all compositions in the considered range. This conclu- 
sion is in agreement with the one made on the basis of 
the inverse temperature hypothesisO discussed in Sec. | 
for alloys with 20, 22 and 24 at.% Pd, but it does not 
allow the change of the temperature behaviour predicted 
by this hypothesis for the Cu-28 at.% Pd alloy. 

The next step is to check this prediction using the ex- 
perimental data discussed in Sec. [II. These data are 



quite limited, since sets of the SRO parameters are avail- 
able only for three compositions. They are also not of suf- 
ficient accuracy for the reproduction of the fine structure 
of the (110) intensity peak. Only in one case, that of the 
alloy 3, the corresponding set is good enough (i.e., con- 
tains sufficient number of the reasonably accurate SRO 
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FIG. 8. Profiles of the off-diagonal part Eod(k) of the AE 
self-energy for the alloys 1 (a), 2 (b), and 3 (c) along the (hlO) 
line. Maximal possible number Na of coordination shells (Ta- 
ble PI) was used in each case. 



parameters) to reproduce the experimentally observed 
splitting in the recalculated diffuse intensity (Fig. ||). 
Even in this case, the recalculation changes noticeably 
the magnitude of the splitting (see Table |). It seems 
that in this particular situation of the split intensity 
peaks even larger sets of the more accurately determined 
SRO parameters are necessary. Nevertheless, the avail- 
able sets can still be used to obtain information about 
the temperature dependence of the peak separation. The 
straightforward approach to this task described in Sec. || 
is applicable only to the alloy 3. The solution of the in- 



7 



I 



in 

I 



I 



21.8 at.% Pd 



(a) 



1073 K 




998 K 



I 



CO 
I 



I 



29.8 at.% Pd 



(c) 




673 K 



> 




> 




k (r.l.u.) 



k (r.l.u.) 



FIG. 9. AE self-energy (a,c) and the related function f{k) defined by Eq. (4.2) (b,d) for the alloys 1 (a,b) and 3 (c,d) 

along the (hlO) line at several equidistant temperatures. Indicated are maximal and minimal temperatures; the temperature 



steps are 25 K (a,b) and 50 K (c,d). The KCM values for the self-energy (Eq. (^) are E^*^^ = -5.87 (a) and E^*^" = -4.78 
(c). 

tive to the accuracy of the {aimn} set than the profile of 
the split intensity peak itself; there is no special reason 
for the self-energy to have any extrema away from the 
special points. Also, Eod is of the second order in aimn 
(Eq. ( ^.Ic )) and therefore decreases in the direct space 
faster than the PCF, which means that the distant SRO 
parameters are less important for its calculation. The 
results of such calculation are presented in Fig. ^. The 
convergence of the results with respect to the number of 
coordination shells included in the AE approximation im- 
proves rapidly with increasing concentration; to achieve 
very good convergence, about 40, 20 and 5 shells are 
necessary for the alloys 1, 2 and 3, respectively. In all 
three cases So( j(k) has a maximum at the (110) position. 



verse diffuse scattering problem given by Eqs. ( |2.12| ) and 
( ^.13 ) cannot be obtained for the alloy 2, because the 
recalculated diffuse intensity becomes negative (Fig. ||); 
in this case the inverse intensity /^^^^(k) and, therefore, 
the effective interatomic interaction VaeO^) would con- 
tain unphysical singularities at those positions in the k- 
space where the diffuse intensity vanishes. The inverse 
problem can be solved for the alloy 1, for which the re- 
calculated intensity is always positive. However, in this 
case the resulting interaction which follows the shape of 
the intensity does not have a split minimum at the (110) 
position. The profiles of VaeC^) ioi the alloys 1 and 3 
are shown in Fig. 0. 

The easiest quantity to calculate in the framework of 
the AE theory of SRO is the direction of the shift. This 
can be done for all three alloys, despite problems with 
the data for two of them as indicated before. Accord- 
ing to Eq. (4.1), the direction of the shift is determined 



by the reciprocal-space behaviour of the self-energy, and 
the l atter can be easily obtained by Fourier-transforming 
Eq. (2.1c), i.e., calculating Sod(k), and using the experi- 
mental values of the SRO parameters. It is expected that 
the off-diagonal part of the self-energy is much less sensi- 



and from Eq. ( O ) it follows that the intensity peaks are 
shifted towards this position (see Fig. This result is 
in agreement with our interpretation of the discrepancy 
between the experimental and the KKR-CPA values of 

q- 

The actual value of the shift can be calculated only 
for the alloy 3, since for the other two alloys positions 
of neither peaks of the recalculated intensity nor min- 
ima of the AE interaction are available. The 10-shell AE 
approximation was used for this and all other calcula- 



8 



o „ 
K LO 



29.8 at.% Pd 



673 K 



0.5 




k (r.l.u.) 



FIG. 10. AE intensity profiles for the alloy 3 along the 
(hlO) line at several equidistant ternperatures. The tempera- 
ture range and step are as in Figs, ^c) and ^(d). 

tions discussed in the rest of the paper. The result is 
0.018 r.l.u.; it should be compared with the deviations of 
the experimental points from the GS line in Fig. ^ The 
deviation for the alloy 3 calculated by the linear inter- 
polation of the GS results is 0.022 r.l.u., which is very 
close to the result obtained in the AE calculation. If we 
assume that the AE shift is about the same for both the 
experimental and recalculated intensities, then the posi- 
tion of the VaeO^) minimum for the former just falls on 
the GS line. The deviations of other experimental points 
are of the same order of magnitude. 

The change of the splitting with temperature can be 
analysed for the two cases (alloys 1 and 3) in which the 
inverse problem can be solved. This is done by calculat- 
ing the self-energy as a function of temperature. Let us 
consider the profile of the self-energy along the (hlO) line. 
Along this line the self-energy is a function of just one 
component k of the wavevector k = (fc, 1, 0). We define 
a function 



(4.2) 



the temperature dependence of which, according to 
Eq. ( |4.1| ) , determines that of the splitting. The functions 
S(/c) and f{k) at different temperatures for the two al- 
loys are displayed in Fig. ^ Accuracy checks show that 
the 10-shell approximation works very well for the alloy 3 
and is still satisfactory (though noticeably worse) in the 
case of the alloy 1. Note that the AE results for the 
self-energy differ considerably from its KCM values; the 
KCM expression for the self-er ierg y can be ob tained, e.g., 
from the comparison of Eqs. ( ^.9| ) and ( ^.14 ): 



^KCM 



c(l - c) 



(4.3) 




In both cases the absolute value of /(fc) decreases with 
increasing temperature for any given value of k, which 



900 



FIG. 11. Intensity peak position q as function of temper- 
ature for the alloy 3 (filled triangles, solid line). Results for 
three other compositions (open symbols, dashed lines), 20 
(circles), 24 (squares) and 28 (triangles) at.% Pd, calculated 
using the same AE interaction are also shown. 



corresponds to the increase of the splitting with temper- 
ature. This result agrees with the conclusion based on the 
assumption of the monotonic temperature dependence of 
the splitting which was made earlier in this Section. 

Finally, a quantitative calculation of the temperature 
dependence of the intensity peak position q can be carried 
out, as before, for only one composition (alloy 3). The 
corresponding intensity profiles are presented in Fig. [lC| , 
while Fig. |l^ shows results for the function q{T) for this 
composition, as well as for three other concentrations 
covering the interval which was considered in Ref. ^. In 
all calculations the same AE interaction (namely, that 
obtained by solving the inverse problem for the alloy 3; 
see Fig. 0(b)) was used. The aim was to find out whether 
the variation of composition for the same interaction 
would lead to any particular change of the function q(T). 
No such change takes place, as it is seen from Fig. |l^; 
the splitting increases monotonically with temperature 
for all four alloy concentrations. 

In summary, we studied theoretically the temperature 
dependence of the Fermi surface-induced splitting of the 
(110) SRO diffuse intensity peak for Cu-Pd alloys under 
equilibrium conditions. The comparison was made with 
experimental observations for these alloys in a steady 
state under irradiation. The validity of the inverse tem- 
perature hypothesis proposed previously to relate the two 
regimes was examined. At equilibrium this hypothesis 
predicts the qualitative change of temperature behaviour 
near the CuaPd composition, namely, the increase of the 
splitting with increasing temperature in the composition 
interval 20 to 24 at.% Pd and its decrease with temper- 
ature as the concentration of Pd increases to 28 at.%. 
Comparing available electron and X-ray scattering data 
with the results of the KKR-CPA electronic-structure 
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calculations, we found that the theoretical approach over- 
estimated the experimental splitting. This disagreement 
was interpreted as the result of the shift of the diffuse 
intensity peaks with respect to the positions of the corre- 
sponding minima of the effective pair interatomic inter- 
action towards the (110) position. An additional assump- 
tion about monotonicity of the temperature dependence 
of the sphtting led to a connection between the direction 
and temperature behaviour of the peak shift. Under this 
assumption the shift towards the (110) position is equiv- 
alent to the increase of the splitting with increasing tem- 
perature. For Cu-Pd alloys this means that the splitting 
increases with temperature for all concentrations in the 
considered compositional range. This conclusion seems 
to be confirmed by the AE calculations, which are, how- 
ever, based on limited experimental data. It agrees with 
the prediction of the inverse temperature hypothesis for 
lower Pd concentrations (20 to 24 at.%) but, contrary to 
this prediction, does not allow any reversal of the tem- 
perature behaviour with increasing concentration of Pd. 
It alsOjContradicts the results of recent computer simula- 
tionaiZI and X-ray measurementsQ for higher Pd content 
(25 and 29.8 at.%, respectively), according to which the 
splitting is (almost) temperature-independent. These re- 
sults are consistent, on the other hand, with the reversal 
scenario. Among possible reasons for this disagreement 
are (i) limited validity of the inverse temperature hypoth- 
esis, (ii) insufficient accuracy and/or size of the available 
sets of the experimental SRO parameters, (iii) pair char- 
acter of the interatomic interactions used in the AE the- 
ory of SRO and (iv) approximate character of the AE 
calculations. Further direct measurements of the split- 
ting as a function of temperature (as in Ref. |^) at the 
discussed range of compositions are necessary to clarify 
the situation. 

The authors would like to thank H. Reichcrt and col- 
laborators for communicating their experimental results 
(Ref. ^ prior to publication. 
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